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We present arguments for the hypothesis that under some conditions, triple correlations of density fluctuations 
in fluids can be detected experimentally by the method of molecular spectroscopy. These correlations manifest 
themselves in the form of the so-called 1.5- (i.e., sesquialteral) scattering. The latter is of most significance 
in the pre-asymptotic vicinity of the critical point and can be registered along certain thermodynamic paths. 
Its presence in the overall scattering pattern is demonstrated by our processing experimental data for the 
depolarization factor. Some consequences of these results are discussed. 
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1. Introduction 

The intensity 7 of light scattered by a one-component fluid drastically increases as the critical point 
is approached fr, 'z']. The physical nature of this critical opalescence phenomenon is well-known, i.e., 
an increase in the magnitude of permittivity (in fact — density) fluctuations and the development of 
long-range correlations between them. The result is that 7 is contributed to by not only single scattering 
effects, but also by those of higher multiplicities. One could expect that 7 provides certain information 
about higher-order correlation functions for the density fluctuations. 

However, a common view is that it is only double scattering j 3- 8] together, probably, with Andreev's 
scattering (due to fluctuations in the distribution function of thermal fluctuations) |9-11] and triple scat- 
tering 1 12] that contribute significantly near the critical point. The pertinent theories are based upon 
quasi-Gaussian statistics for the fluctuations. As a consequence, any information on the irreducible parts 
of higher-order correlation functions is lost and multiple scattering is viewed as parasitic. 

By contrast, we pursue the idea that the so-called 1.5- (sesquialteral) scattering, caused by triple cor- 
relations of density fluctuations, contributes significantly to 7. In this report, we support this statement 
by our results of processing extensive experimental data flS*] on the depolarization factor near the criti- 
cal point. In our view, these results strongly indicate that the 1.5-scattering is noticeably present, under 
certain conditions, in the overall scattering pattern. Moreover, the idea of 1.5-scattering allows us to give 
simple explanations for the anomalies in the behavior of 7"^ in the gravity field fli.'lS'j and those of the 
Landau-Plazcek ratio near the A-line 1 16-18], which were observed long ago, but have been interpreted 
controversially. 

It is also important to note that, according to our own theory and in view of earlier estimates f3], 
from among different three-point configurations of density fluctuations 8p (located at points r, ri, and 
r'), the major contribution to the integrated 1.5-scattering intensity 7i.5 is made by those in which two 
successively scattering fluctuations merge on the macroscopic scale (|r-ril 0; technically, this is a 
consequence of the replacement of the internal electromagnetic field propagators in the iterative series 
for the scattered field with their leading short-range singularities). As a result, 7i 5 is determined by the 
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Fourier transform of the pair correlation function ([5p(r)]^5p(r')>. Provided Polyakov's hypothesis fli^l 
(see also |20, 21]) of conformal symmetry of critical fluctuations is valid, the latter is expected to vanish 
at the critical point (see appendix A). It follows that the recovery of the 1.5-scattering contribution from 
/ and a scrutinized study of its behavior along appropriate thermodynamic paths ending up at the crit- 
ical point provide a unique opportunity for experimental verification of Polyakov's hypothesis Hai for 
systems with scalar order parameters. 



2. Theory of 1.5-scattering 

2.1. General expression 

The theory of 1.5-scattering was propo sed in |22] and developed further in (i^ll^l; some additional 
numerical estimates were made in |25l|2^. We assume that molecular light scattering from condensed 
matter is a result of re-emission of light not only by single fluctuations, but also by compact groups of 
fluctuations. By compact we understand any group of fluctuations all the distances between which are 
much shorter than the wavelength A of the probing light in the medium. Physically, scattering by such a 
group is single. The overall polarized single-scattering spectrum is, therefore, given by the series (iill 

/(q,w)= Inm{q,n), (2.1) 

n,m^l 

where 

- — J x-Re j drj dr([<5e(r,r)]"[5e(0,0)]'")e'"^-"! '' (2.2) 

V 

is the contribution from a pair of compact groups of n and m permittivity fluctuations [attributed further 
to density fluctuations, Se =; [d£/dp)jdp], eq is the equilibrium value of the permittivity, Q and q are the 
changes in the light frequency and wavevector due to scattering, and the scattering volume V is included 
into the proportionality coefficient. 

It is only the term In (q, ^) in equation <2.1) that has been associated so far with the single scattering. 
The 1.5-scattering intensity is defined as 7i.5(q, n) - /i2(q,n) -i- hi (q, f^)- 



2.2. Hydrodynamic region, qvc « 1 

Far enough from the critical point, where the correlation radius rc « A and nonlocal correlations 
between fluctuations can be ignored, the integrated 1.5-scattering intensity can be expressed in terms of 
the third moment of thermodynamic density fluctuations Ap | 22i|: 



2 

3co 



de 



V 



T,V 



(2.3) 



where /^r is the isothermal compressibility of the fluid and y is a macroscopic volume over which the 
fluctuations Sp are averaged to single out their thermodynamic parts Ap. We suggest that V is slightly 
dependent on temperature far away from the critical point, but oc oc in the critical region. 
Calculations with the van der Waals and Dieterici equations of states give the estimates 



/1.50c p— ^^-ePcWjfil 

£q\ dpi V 



and 



/1.5 oc 



2 

3£o 



de 



3 ]^2j2 



respectively, where w = pc/p - 1, Iwl « 1 is the deviation of p from the critical value pc and Pc is the 
critical pressure. It follows that the 1.5-scattering can become of significance in those domains in the 
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(T,a»)-plane where co^O, but Pt is sufficiently large. Then, I\,s oc toPj for a non-critical isochore, but 
/i.5 cx jS^^ or even 7i.5 for the critical one. 

A distinctive feature of the l.S-scattering contribution is that it is not positive definite: for instance, 
/i.5 < in the region where ld>0 and t > 0, at least. 



2.3. Fluctuation region, qvc » 1 

Understanding, in this section, p as a scalar order parameter, we see that formulas and <2.2> agree 
with the hypothesis of algebra of fluctuating quantities (i^. Then, within the first order of c-expansion 
and in the long-wave limit q ^ 0, the critical index of 7i.5, defined by 7i.5 oc \t\~^'", is estimated to be 
fi2i ~ 0.67 for 0) - 1 22]. This value is close to an earlier estimate of 0.7 given in |20J. Correspondingly, 
7i.5 oc on the critical isochore and in the immediate vicinity of the critical point. 

This result can be refined using the algebra of fluctuating quantities (see appendix A). However, it 
is more important to emphasize that it disregards the conformal invariance hypothesis |19]. If the latter 
is indeed valid, then the orthogonality relation holds for fluctuating quantities with different scaling 
dimensions (see H^lill]), that is, 7i.5 ^ as the critical point is approached. 



2.4. Intermediate region, qvc < 1 

This region is of special interest to us because it is tjrpical of actual experiment. Taking into account 
that correlations between fluctuations 6p remain relatively weak, we argue (H] that the convolution-type 
approximation 

2 c' 

(PkiPkiPkj)--- 7=G(A;i)G(A;2)G(fc3)5k3 -ki-k2 

can be used for the three-point correlation function of density fluctuations. Here, pk is the Fourier compo- 
nent of 5p(r), G{k) = (|pkl^>. and c' is a k-independent function of temperature and density. Calculations 
with the Ornstein-Zernike expression for G{k) then give: 



c' / deY p'ki.T'-pt iqrc\ 



D^klT^Pl 

Requiring that in the limit qr^ « 1, equation <2.4> transforms into equation d2.3> . we recover c' through 
the third thermodynamic moment of density fluctuations, with an accuracy to a positive proportionality 
constant: 

- 'If )^ J 

Extrapolation of formulas I |2.4> and 1 12. 5> on the fluctuation region shows that (see appendix B) c' ^ 
and, therefore, 7i 5 ^ as both t ^ and &) ^ 0, which is in accordance with the conformal invariance 
hypothesis. 

The structure of the l.S-scattering spectrum in the intermediate region is discussed in i24il . 



3. Depolarization factor 

3.1. Theoretical considerations 

Now, we are in a position to scrutinize the effect of l.S-scattering on the depolarization factor A as 
a function of temperature (in fact, (ij) and the geometrical size L ~ V^'^ (volume V) of the scattering 
system. Suppose that the following contributions to 7 are present in the intermediate region qr < 1: 

(1) the "standard" intensity In oc VPr of polarized single scattering due to density fluctuations (ll01; 

(2) the intensity 7i.5 of polarized l.S-scattering |22-24]; (3) the intensity 7a oc V of depolarized single 
scattering due to anisotropy fluctuations 1 1] (it is virtually insensitive to the critical point); (4) the intensity 
hp oc of polarized double scattering due to density fluctuations (5) the intensity hd oc 
V^'^P^ of depolarized double scattering due to density fluctuations |3-8]; (6) the intensity 7a oc VPj^ 



13003-3 



M.Ya. Sushko 



of depolarized single scattering due to fluctuations in the distribution function of thern\al fluctuations 
(Andreev's scattering) L9J. Then, A is given by 



4 + hd + Ia 

(3.1) 



h + h.5 + hp 



In view of the individual temperature dependences of the above contributions and under the condi- 
tion /i.5 = 0, A as a function of fir is expected to decrease first, then reach a minimum, and then increase 
again. Such a behavior, indeed observed in the experiment, is considered as a manifestation of double 
scattering effects. However, as we show later, the presence of the 1.5-scattering contribution does not 
alter this qualitative behavior of A as a function of fir- 

Thus, expression d3.1> should be transformed in order to obtain an experimentally-measurable func- 
tion whose behavior significantly depends on whether the 1.5-scattering contributes to A or not (iill . 
Rewriting <3.1) as 

_ i+/i.5/r^+-f2p/r^ 



(3.2) 



and taking into account the specific features of the intensity contributions, we immediately arrive at the 
relation 



A l + cL-^p-^ + dL-^p- 



(3.3) 



valid for the intermediate region qr <l. Here, the coefficients b, c, and d are practically temperature- 
independent and positive constants. The coefficient a oc c'arctan(^) /rc is due to the 1.5-scattering 
contribution and is not positive definite. If the 1.5-scattering is negligible, then a-0 and the right-hand 
side in formula B.3t is a monotonous increasing function of fij- With the 1.5-scattering present, this 
monotonous behavior is expected to be violated. The effect should be most pronounced in the following 
two cases. 

(1) The critical point is approached along a noncritical isochore w > 0. Then, /i 5 oc -j6^'^ and a is 
close to a negative constant. 

(2) The critical point is approached along the path where t ^ 0, w ^ 0, and 7i.5 > 0. The relative 
magnitude h.slh of the 1.5-scattering should start decreasing somewhere due to the temperature law 
/1.5 oc P^j^ coming into play 1 20, 22] (see section l23) or as a consequence of the conformal invariance (3 
[2^. As such a path, the liquid branch of the coexistence curve can be quoted. 

Thus, by varying the temperature ifij) and density (oi) of the scattering system, we hope to "stick out" 
the 1.5-scattering contribution from among the others. It should manifest itself as a non-monotonous 
behavior of the experimentally-measurable quantity L/^j-ZA"^ with fij- The fact that the scattering con- 
tributions involved depend differently on L, provides an additional powerful option for analysis. Some 
results obtained by processing the extensive depolarization factor data (Hi are presented in figures[TtiT2] 
They generalize our earher results Esll . 



3.2. Data processing 

3.2.1. Noncritical isochores p < pc.r >0 

Figure [1] represents the log-linear plots of the quantity L(_DA)~^ oc LjSj-A"^ as a function of the pa- 
rameter (fco is the wavevector in vacuum of the incident light) 

D-^ = (eo - 1)^ (£0 + 2f kB TPt oc Pt 

for xenon along the a) - 6.8 x 10"^ isochore and five values of L (in cm). The parameter D (in m) is 
a convenient measure of the distance to the critical point |27j|. It is evaluated in |13] as a function of 
temperature and density by using the Clausius-Mossotti relation for eq and the restricted linear model 
equation of state |28] for Pt- These calculations are claimed to be most reliable for the region not very 
close to and not far away from the critical point, i.e., the one of special interest to us. 
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Figure 1. L(DA) ^ versus D ^ along the w = 6.8 x lo ^ isochore of xenon for five values of L, based on 
experimental data . From left to right, r decreases from 7.8 x lo ^ to 1.2 x lo ^ . Three segments can 
be distinguished on each of these log-linear plots. 

As T reduces from 7.8 x 10~^ to 1.2 x 10~^, three typical temperature imervals are clearly seen in 
figure[T] We shall refer to them as segments A (leftmost), B (central) and C (rightmost). It is easy to note 
that simple division of L(Z)A)~^ by L and changing to (DA)~^ as a function of for different values of 
L transforms the original plots dissimilarly on these segments: the plots tend to merge on A and C, but 
invert the vertical ordering and disperse on B (figure|2). This fact imphes that (DA)~^ is contributed to 
by terms with differing functional dependences on L. Our further analysis of them is guided by relation 

Suppose that on segments A, i.e., the most distant from the critical point, 7i prevails much over /i s 
and l2p. Then, relation 1 13. 3> takes the form 



A 



1 



l + cL-^p-^ + dL-^P 



-3/2 ■ 



It follows that the dependence of (DA) upon should be close to a linear one, with the slope independent 
of L and, if the Andreev contribution is noticeable, a slight concavity: (DA) oc const + cD^ + d[D^)^^'^. 
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Figure 2. (DA) versus D for the data shown in Figure 3. (DA) versus D for segments A. From right 
figure[l] to left, t decreases from 7.8 x io~^ to 3.3 x io~^ . 
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Figure 4. L(DA) ^ versus D for segments A and 
B. From left to right, t decreases from 7.8 x 10~^ to 
1.4 X 10"*. 



Figure 5. L(DA) ^ versus D ^'^ for the portions of 
segments B where t decreases from 9.8 x 10~* to 
1.4 X 10"*. 



Figure[3]does exhibit, at least approximately, such kind of dependence. The study of the latter could, in 
principle, provide experimental estimates for the magnitude of 7a- However, the discussion of this issue 
is beyond the scope of the present report. 

On segments B, where we expect l2d to dominate over /g and 7a. but Izp to remain relatively weak as 
compared to h and /i s, relation takes the form 

A 

At w > 0, L(DA)"^ should decrease with D~^'^ by the linear law L(7)A)"^ oc -D"^'^, with negative and 
equal slopes for different values of L (figures|4]and[5). Correspondingly, the dependences of (DA)"^ upon 
-73"^'^ should be linear on B, with slopes inversely proportional to L, but, as was already mentioned, 
merge on A (figure [6). 

Finally, segments C are formed mainly by single and true double scatterings. Relation <3.3) takes the 
form 

LBt 

oc 1 + bipr 
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Figure 6. (DA) ^ versus D ^'^ for the data shown in figure|4] 



13003-6 



Observation of triple correlations in fluids 



700 r 



600 
500 
400 



300 



L, cm: T 0.547 

□ 0.313 

+ 0.169 

o 0.035 

▲ 0.021 



100 [ail ^ 



(DA)-' 



2500 



1000 ) 



L, cm: T 0.547 

□ 0.313 

+ 0.169 

o 0.035 

A 0.021 



ko4iV tat, 



500 600 



100 200 300 400 500 600 



Figure 7. L(DA)~^ versus for tfie entire w = 
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Figure 8. (DA) ^ versus D ^ for the entire w - 
10~^ data plotted with a linear scale. 



6.8 X 



and we expect L(DA)"^ to increase with by a linear law, with slope proportional to L (figure [7). 
The (DA)~^ versus plots for different values of L should approach a single straight segment as 
increases (figure |8). 



3.2.2. Liquid branch of the coexistence curve 

The dependence of L(DA)~^ upon along the liquid branch of the coexistence curve of xenon is 
shown in figure[9l It agrees well with our expectations. 

Thus, the above processing of experimental data (Hi clearly reveals the presence in the overall scat- 
tering pattern of a contribution which we associate with the 1.5- (sesquialteral) molecular light scattering. 
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Figure 9. L(DA) ^ versus D ^ along the liquid branch of the coexistence curve of xenon; L = 0.547 cm. 
From left to right, t changes from -9.2 x 10~^ to -2.9 x 10~^. 
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3.3. Numerical estimates 

Now, we present quantitative estimates of the magnitude of l.S-scattering intensity. They were ob- 
tained by fitting the L(DA)~^ versus data for the entire w = 6.8 x 10"^ isochore and then used to 
reproduce the original A versus D data fl^ . 

In view of formulas | |3.2> and <3.3> . the fitting function was taken in the form 

fix)^ ^ [K+Ax^'^ + Bx), x = D~\ 

and the foUow/ing two sets of coefficients were chosen: Ki = 581.066, Ai = -63.6968, By = 3.66569, d = 
0.012 and K2 = 607.508, A2 = -89.0932, B2 = 5.79532, Cz = 0.012. The relative magnitudes ri.5 = h.slh 
and r2p = hplh of the l.S-scattering and polarized double scattering, as compared to the single scattering, 
were estimated as 

Ax^''^ Bx 
ri.5 = ^, r2p = — . 

To calculate A with //, our theoretical estimate (llll hd ~ \hp was additionally used. Then, 

Bx 

A= . 

&Kf{x) 




Figure 10. Fitting the (DA)~^ versus data (t) = Figure 11. Fitting the A versus D data (w = 6.8x lO"-*, 
6.8 xio~^,L = 0.547cm) with /i (dashed line) and /2 L - 0.547cm) with /i (dashed line) and /2 (solid 
(solid line). line). 

The results are demonstrated by figures [T0lfT2] They clearly show that in the intermediate region, 
the intensities /i 5 and Tap reach magnitudes comparable with that of I\, but are opposite in sign and 
tend to compensate for each other. These facts are surprising. They contradict the common view that 
multiple scattering contributions come into play gradually as the critical point is approached. In other 
words, they imply an asymptotic nature of the iterative series for the overall scattering intensity near the 
critical point. They can also be interpreted in the sense that triple and quadruple correlations in fluids 
contribute, at least to light scattering effects, in opposite directions. 

The agreement of our fitting results with the A-data (l^ (figure [TT) is also unexpectedly good. 



4. Conclusion 



The above estimates have prompted us to identify other experiments where the situation is favorable 
for l.S-scattering to come into play. First of aU, of interest are the studies on light scattering from critical 
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fluids under the earth's gravity. Due to the gravity effect, the system is spatially inhomogeneous in the 
vertical direction. A negative 1.5-contribution is expected to appear in the hght scattered from the fluid 
layers located above the level of critical density. For such a layer, the total scattering intensity / should, as 
a function of temperature, start decreasing somewhere as the critical point is approached. In other words, 
the versus t dependence should have a minimum at some tq. The effect was indeed registered, for 
instance, in freon 113 jl3,Ell- Our estimations for the minimum location, tq ~ 10"^ for heights up to 
20 mm (t > 0), agree well with experiment. 

We suggest that a specially-designed processing of the gravity-induced height- and temperature de- 
pendences of 7, obtained for systems with a scalar order parameter, is a feasible opportunity for singling 
out the l.S-scattering contribution and verifying Polyakov's conformal invariance hypothesis. 

To finish, we mention that the studies of the spectral distribution in critical opalescence spectra are 
of great interest as well. In particular, we have proved that in the presence of l.S- and double scattering 
effects, the ratio of the integrated intensities of the Rayleigh and Brillouin components takes the form 



, 1 + ai.5 ri,5 + (22p'"2p 
l + ^'l.5'"l.5 + fo2p'"2p ' 



iJexp^i?^^^^^^^^^^, (4.1) 



where i? = y - 1 is the well-known Landau-Placzek ratio (sofl for single scattering (j = cp/cy, cp and cy 
being the specific heats at constant pressure and volumes). The coefficients ^1.5 and are given in fi^l . 
whereas fl2p and b2p can be recovered from the results l3lll : 

a2p^l-^+—^, fo2p = 2-^. (4.2) 
27 2(7-1) ^ 2j 

Suggesting that ri.5 - 0, it is not difficult to verify, based on experimental data fii*] for He^, that the double 
scattering alone should cause i?exp to exceed R as the A-line is approached along a high-pressure isobar. 
Such a fact was indeed registered (l8|], but in most other experiments in this series the tendency was 
direct opposite |16-18]. We attribute the reduction in _Rexp to the effect of l.S-scattering. 
Our detailed calculations of the above effects will be presented elsewhere. 
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Appendix A 

Let Afc(r) be a complete set (algebra) of fluctuating scalar quantities with scaling dimensions Aj;: un- 
der the scaling transformations r Ar, AkiAr) A~^* Afc(r). According to the local algebra hypothesis 
(see, for instance, i2oll ), the scalar function 5p(r) and, therefore, the scalar function [(5p(r)]^ can be de- 
veloped into the series 

oo oo 

5p(r) = ^ flfc^fcCr), [Sp{r)] = ^ bmnBmnir], 

k=l m,n=l 

where the coefficients bmn are simply related to the coefficients ajc, and BmnM - ^m(r)^«W are scalar 
quantities with definite scaling dimensions Am„ - Am + A„: Bmni^r) A''^'"" BmnW- Correspondingly, 
the pair correlation function {Sp(r'] [5p(r)]^> is given by a linear combination of pair correlators 
Kk,mnQr:'-r\) = {Ak[r')Brr,„{r)). 

For a rf-dimensional system with scalar order parameter, only the Fourier transforms of those Kk,mn 
can reveal a singular behavior near the critical point which satisfy the condition A^-i-Amn < d. In the first- 
order c-expansion Afc = jk{2-e) + ^k{k-l)£, where c = 4 - d |20]. If = 3, then the relevant correlators 
are ^Ti^i, ^^1,12 = ^1,12, and Kz^n, each involving two scalar quantities with different scaling dimensions. 
Once Polyakov's conformal symmetry hypothesis holds true and the system is spatially homogeneous 
and isotropic, they all vanish at the critical point due to the orthogonality relation (see (2^ [2l|] and the 
literature cited therein). 



Appendix B 



As one of the ways for evaluating c' in the fluctuation region, we can use the asymptotic equation of 
state 

which immediately follows from the requirements that it leads to (a) the correct asymptotic behavior 
of a limited number of the fluid parameters along the selected thermodynamic paths (the susceptibility 
%- P^Pt along the critical isochore, the critical isotherm equation, and the derivative of pressure with 
respect to temperature at the critical point) and (b) reveal a Van-der-Waals-type loop below the critical 
point. In this equation, n - P/Pc -1,t- T/Tc-1,ci)- V/Vc - 1, To and Do are the critical amplitudes for 
X and the critical isotherm, respectively, and M is a constant. The definition 

v[dP)T.v Pc[l + co)[dn 

relations 

'dPr] (dPr] (dPr] (dV] (dPr 



oP It,v I dP It,n I oV lT,N\dP It.n \ 00 



and formula <2.5> then yield 



-PcDoS[S-l)p~^Ci)\wf~^. 
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EKcnepi/iMeHTaiibHe cnocTepe>KeHH5i norpmHi/ix Kopeyini^m 
y piAi/iHax 

M.51. CyuJKO 

OpfiCbmvi HaMioHa/ibHMM yHiBepcMTeT iMeHi I.I. MeMHUKOBa, By;i. /iBop^HCbKa, 2, 65026 O^eca, YKpaiHa 

HaBefleHO apryMeHTH Ha KopucTb rinoTesi/i, no npn neBHUx yMOBax MeTOflOM Mo;ieKy;inpHoT cneKTpocKoniT 
MOXHa sapeecTpyBaTM noTpiiiHi Kope;ifli4iT <j);iyKTyai4iM rycTMHH b piflUHax. I^i Kope;ini4iT npoflsyiflKSTb ce6e y 
BMryiflfli TaK ssaHoro 1.5- (to6to nisTopaKpaTHoro) poscinHHfl, s\Ke e HaM6mbLU cyrreBi/iM b nepeflacuMnroTi/m- 
h\\a o6;iacTi KpnTW4H0T tohkm ra Moxe 6yTW sapeecrpoBaHe BSfloex nesHnx TepMOflWHaiviiHHUx LU/iaxiB. Moro 
npucyiHicTb y sarayibHiw KapinHi posciaHHH fleMOHCTpycTbcn peayjibTaraMi/i o6po6KM b\p,om\ax eKcnepuMen- 
TayibHMx AaHMx p^!\ Koe(t)iLiieHTa Aeno;inpM3aL4iT. 06roBopeHo fs,eaK\ Haa\\fl,m uyix pe3y;ibTaTiB. 

lOiiOMOBi cnoBa: tpjiyKTyaiii'f rycrnHH, KprnmHa ona/ieaieHLiin, poscinHHH KpsTHOcri 1.5, KoetpiifieHT 
ftenompiAsaiiii 
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